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' Abstract 

in : 

We compute a compact OPE formula describing power corrections to the perturbative 

■ — % — 

f^) ■ expression for the asymmetric MOM-renormalized running coupling constant up to the lead- 

| ing logarithm. By the use of the phenomenological hypothesis leading to the factorization 

of the condensates through a perturbative vacuum insertion, the only relevant condensate in 
I \ the game is (A 2 ). The validity of the OPE formula is tested by searching for a good-quality 

C^H 

| coherent description of previous lattice evaluations of MOM-renormalized gluon propagator 

and running coupling. 
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That the running coupling constant can be extracted from the three-gluon vertex in the 
Landau gauge was proposed several years ago in a seminal work Q. This method appears closer 
to a physical interpretation than Scrodinger functional's ^ but, and mainly for the same reason, 
more systematic effects should be managed to produce a reliable prediction for the perturbative 
as- In particular, the first statistically meaningful attempts to follow that method missed the 
impact of power corrections 0] and failed to give an estimate for the coupling constant comparable 
to others in literature (that of ref. Q, for instance). In parallel, these power corrections to the 
gluon propagator and to the running coupling constant have been largely studied in the last 
years gg. 

Although the first trials were rather inconclusive ||, momentum power contributions have 
been manifestly put in evidence [J5],||] for the Landau gauge three-gluon coupling constant renor- 
malized in both symmetric (MOM) and asymmetric (MOM) momentum substraction schemes. 
In ref. M the parameter A^jg is estimated from the matching of MOM as lattice results to 
a perturbative three-loop formula corrected by an unavoidable l/p 2 -term. This naive ansatz 
used in ref. j/J seems to eliminate most of the systematic deviation from the three-gluon vertex 
estimate of the perturbative as and a precise prediction of A^jg emerges in full agreement with 
that of ref. Q. Unfortunately, the errors quoted in this work were clearly underestimated. In 
fact, the prediction of A^? is so sensitive to a logarithmic dependence on the momentum scale 
of the coefficient of 1/p 2 that it appears to range over an interval of 40 MeVg. 
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On the other hand, in refs. [||,|9) a description in terms of OPE for the power corrections to 
MOM three-gluon as is successfully tried consistently with an analogous description for gluon 
propagator. The OPE approach provides through SVZ factorization [ 10 1 a perturbative tool to 
obtain the leading logarithmic dependence of non-trivial Wilson coefficients. In particular, it 
can be applied to compute the coefficient of A 2 in the MOM case for Landau gauge three-gluon 
coupling constant ||. Obviously, the expectation value of A 2 in a non-perturbative vacuum is 
also estimated via this OPE approach. A non gauge invariant condensate, that, being in Landau 
gauge, is invariant under infinitesimal gauge transformations and thus connected in some way 
to the gauge invariant (-A^in) defined in [11], is then computed. 

The same approach for MOM coupling constant requires to use 



T* A a J-p)Ai(p) 



' fj.ua' • 



+ (P) :<'(0M fe ;,(0)<,(0): + ... ; (1) 



where the expansion to the three-gluon local operator is necessary because the three-point Green 
function in MOM can be written as 



( T* (2»(-p)It(p)I°(0)) )np ss 

(0|T* (A°(-p)A b v (pj) |<# NP = ( Cl )^ (0| : A$(0) : |^)np 

+ (0|:<,(0)^',(0)^;(0):| 5 p NP + ••• .(2) 

The index NP refers to the non-perturbative nature of the vacuum state in Eq. (|2|), while T* 
refers to the standard time ordered product in momentum space. It should be noticed that: i) 
No other local operators with the same dimension of A 2 are written in Eq. (|l]) because, unlike 
the identity or A 2 itself, they do not generate non-null vacuum expectation value (v. e.u.)[] ii) 
Operators other than the three-gluon local one for the same dimension, as d^A^(0)Ap(0), do 
not appear explicitly because they are phenomenologically supposed not to survive, as will be 
argued below. The identity and A 2 clearly do not contribute to the matrix element considered 
in Eq. (D 

Following standard SVZ techniques to obtain the peturbative expansion of OPE Wilson 
coefficients, we compute the appropriated matrix element of Eq. (jl|)'s l.h.s. to the wanted 
order. It is immediate that taking the perturbative vacuum leads to 

<0|T* (l»(-p)j5*(p)j5J(0)) |0) = {pCi% . (3) 

where c\ can be straightforwardly identified with the perturbative expansion for Green Function 
with an amputated soft-gluon leg, 

rg£(p,-P,Q)^ -2^V»/ a a' G%(p 2 ) , (4) 

x No Lorcntz invariant tensor with an odd number of indices can be built without non-zero momenta. 







which should be proportional to its Landau gauge tree-level tensor. Beyond this purely pertur- 
bative first contribution, the situation is a bit more complicated. We consider now the matrix 
element in the l.h.s. of Eq. (||) between two external soft gluons and the perturbative vacuum. 
Then, we obtain up to the considered order in perturbation theory 



(0\T* (A°(-p)Al(p)A C p (0)) bUUnected 

= <0| : A$(0)J& (0)4(0) : \g c p g[g s a ) 



i pva 

For the matrix element in the r.h.s. we have 



(5) 



(0|:^(0)4'(0)4(0):|^i^) 

• ss' W ct f 1 

= G(2) CTCT ,(0) G(2) AV (0) G<V(0) [V s Xl gj'g*g p n a ' sl 5 b ' v 5 c, t + O(a) } , (6) 

where V refers to all the possible permutation of the couples (cr's'), (A'/'), (rt). The Wilson 
coefficient C3 may thus be computed at tree- level order as: 



Kil (c 3 ) 



_ (A[(0)A-(-p)Al(p)A c p (0)AU0)) 

pus'l't ss' IV ct 

G(2) ff(T ,(0)G(2) AV (0)G(2) pr (0) 



(7) 



The ratio in Eq. (|7|) represents symbolically all the tree-level diagrams with five gluon legs where 
the three of them carrying zero momentum are cut (See fig. fi]). 





Figure 1: Diagrams contributing to the tree-level Wilson coefficient in Eq. (Q). Crosses mark the soft 
gluon legs coming from the condensate. 

Had we directly dealt with the three-gluon local operator : A p ,(0)A b ^,(0)A c p ,(0) :, the evalu- 
ation of the higher order tensors involved would require a much more tedious calculation. On 
the other hand, the vacuum insertion between one gluon field and the other two leads to 



with 



(0[:AV0)j£(0)A£(0):|rf) 



NP 



4(JV? _ i) 



q-a'b'c'r 
1 p'v' p't 



(8) 



Ta'b'c'r _ ™' b' r re' i „ ra' c' v' r b' , „ XC'b' „ rra' 

V^'p'i V^'° 9p>o t + g^'p'O g T d t + g p > u >d g p ,d t . 



za'b 



c'b> 



(9) 
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We may phenomenologically assume vacuum insertion to work for a certain renormalization 
momentum scale. Thus, the following replacement 



• A 2 At (f)\ ■ 

:A% (0)4', (0)^,(0): - ' ^ ^ (10) 

can be done for Eq. (^), the other components of the tensor in l.h.s. not being required for our 
purposes. Furthermore the local gluon field -A^(O) is to be contracted with the zero- momentum 
gluon field defined for the vertex. The same vacuum insertion assumption leads us to argue 
that, for instance, 



(0|: ^4', (0)^(0) :\g c p ) 



(d'X'WW <%(0,M 



. 



(11) 



Then, the only non-zero surviving condensate comes from the three-gluon local operator. 

It is easy to see that, using Eq. (|io|), the relevant coefficient multiplying the local operator 
in Eq. © is 



(r) a b ^ P ' (n)ra ^ T 1 K(0K(-p)4(pK(0)^(Q)) v 

am a ao) gw Xv , (o)gw_(o) 



where the r.h.s. may be straightforwardly obtained from Eq. (|7]). This last Eq. (12) gives 
the prescription in which Lorentz and color indices of external gluon legs are contracted in the 
diagrams contributing to the tree-level Wilson coefficient (see Fig. [l]). 

Thus, since Eq. (H)'s is the only Landau gauge tensor for the asymmetric three-gluon ver- 
tex]^], we can write for p 2 = —k 2 



k\T* A a J-p)Al(p)A c J0)) ) NP = -IfgiMf 



x ^ (^j G(% p (0) + c 3 (g,^) j . (13) 

We know the scalar coefficients ci,C3 in Eq. di~3| ) to be dimensionless by OPE power counting 
rules, and hence both only depend on the bare coupling g and the dimensionless ratio of momen- 
tum over regularization scale Q The term c\ can be obtained from Eq. (|!|) and, at tree-level, C3 
is to be computed by projecting over the Landau gauge tree-level tensor in Eq. (Q) the result 
fromEq. @, 



(^ 2 \ _ / \abfj,'u'p'/ \ rra'b'dr 1 _L iiv t \t t 

= 3 5 . (14) 

If we renormalize following the MOM prescription at momentum scale /i 2 , and then apply the 
assumed vacuum insertion factorization, we can write: 

2 The dependence on regularization momentum scale (a -1 in lattice regularization or e _1 /x in dimensional, for 
instance) has been up to now omitted to simplify the notation 



A 



( T* [A a (—p)A b u (p)A c {Q)) ) 



ci 



A; 2 



R,p 2 



r,a(/x) 



C3 



,a(/x) 



<^ 2 >rV 



, tc 

* G 2 



(0,M 2 ) 



4(iV 2 - 1) -k 2 



(15) 



where, after renormalization, Wilson coefficients depend on the ratio of the momentum and 
the renormalization scale, and on a(fi) = g\(n) / '(4vr) , i.e. the coupling constant consistently 
renormalizecHn^MOM. 

In the MOM scheme, renormalized Green functions take formally, at the renormalization 
scale, the same tree-level value but in terms of the renormalized coupling constant instead of 
the bare one. This defines Z MOM (fi) = fi 2 G^ 2 \fj?) to renormalize appropriately the two-point 
Green function, where is the scalar factor of the bare two-point Green function defined in 
ref. ||. The three- gluon Green function is then renormalized dividing by (Z MOM (/x)) 3 / 2 . One 
gets 



(3), 



g R (k 2 



k 4 ° R 



c 2 V) 



(tfGg>(*V)) 



-1/2 



(16) 



The renormalized scalar factor for the three- gluon Green function with an amputated soft-gluon, 
Gft /Gr , can be projected out from Eq. (15), similarly to what is done for C3 in Eq. (0), 



while for G 



(^ 2 ) 



we can write 



k> Gg>(*V) = c g.aw) + c 2 (^,a(M)) ^ , (17) 

where the scalar coefficients Co and C2 can be derived from those in Eq. (||) and computed 
similarly to c\ and C3, as explained in ref Q. Thus, after replacing in Eq. ( |i~6| ) we will finally 
get: 



gn(k 2 ) = ci 



1 + 



k 2 



-1/2 



1 ' i 



Rj A 1 



1 



4(iV 2 - 1) -A: 2 



(18) 



The purpose is now to compute to leading logarithms the subleading Wilson coefficients in 
Eq. (|i~5|), as done in ref. ||. It will be, to this goal, useful to consider the following operator 
expansion, 



5=5 Pr^ipiU* T* (A^-p)Al(p)A c p (0) 



9k 2 



R lA t 2 



C3 ( 771, 



:^ 2 4(0):^(0) 



R,/i 



+ 



(19) 



where dots refer to terms with powers of 1/k other than 6. If the vacuum expectation in the 
non-perturbative vacuum is considered for the r.h.s. of Eq. (|l9|), the result under the vacuum 
insertion hypothesis is known to be diagonal in color and Lorentz spaces. Then we will contract 
in both sides of Eq. (IS) with g pT 5 c t and we will take the following matrix element (see ref. [112]) 



W P p 9 ± " / (p)fabc (0\T* (Al{-p)Al(p)A%0)) Is^iW 



c 3 (^,a0")) 



(Q\:A'AP(0)--\9o9 s A)w + ••• I (20) 



-k 6 



where we take external gluons in a perturbative vacuum and carrying soft momenta. From Eq. 
© we get 



2 4 P p g^(p)fa b c {0\T* (A»(-p)^(p)A«(0)) \g s J x ) 
~9 k (0| : A*aW) : \gcgsj x ) 

= Z 3 (» 2 )Z-W)c 3 (^,a(»)\ = Z- l/ \^)Z- 1 C3 ^,a(/x)) (21) 
where, always in MOM prescription, Ar = Z% A and Z43 is defined such that 



: A*AP(0) : A c p (0) 

while Z = 3 ^ 2 is an useful notation for the constant renormalizing the matrix element for 

the three-gluon local operator where the external soft gluons are explicitly cut. If we recover 
the divergent factor Z = Z^Z^ introduced in ref. || for the matrix element of two-gluon local 
operator coming from proper vertex corrections, Z can be thought to be decomposed as 



2 = Z- 3 V)^ V V) :^(0):^(0) ; (22) 



Z(n 2 ) = Z(^ 2 ) Z K (fi 2 ) , (23) 

where Z takes the divergent part coming from the diagrams for the matrix element in r.h.s. of 
Eq. ( p0[ ) which can be factorized (diagrams (a,b) of fig. 0) as in Eq. (||), and Z K should be 
computed from those which can not (diagrams (c,d) of fig. |2[). 

Then, taking the logarithmic derivatives with respect to fi in both sides of Eq. (]2l|), we get 
the following renormalization group differential equation: 



-27 (a(ji)) - 7 (a(jj)) + q^j. + p ( a ( M )) A I C3 f^,a(fi) ) = ; (24) 



with the formal solution (see ref. 



/i 2 ' / • V (, (/') 





Figure 2: All the possible leg permutations from diagrams in the figure contribute to the renormalization 
constant Z defined in the text, (a) and (b)-like diagrams, which do not break the assumed-to-work factor- 
ization hypothesis, are renormalized by Z previously computed in Q. (c) and (d) breaking factorization 
diagrams give Z K . The local operators are drawn as gray bullets, the two of them joined by a dashed line 
represent the ones contracted to give A 2 in Eqs. through the replacement in Eq. @. 



where 





d 




dliifj, 2 


7(a(/i)) = 


d 


lnZ 3 (/i 2 ) 


din fj, 2 


/?(«(/*)) = 


d 


a(n) = 



a(/i) 



7o- 



4ir 



+ 7i 



4ir 



+ 72 



4ir 



+ 



(26) 



The boundary condition of Eq. (24) is given by our MOM- like prescription for the renormaliza- 
tion of the condensate by Z43 : the condensate is renormalized such that the Wilson coefficient 
takes the tree-level form at the renormalization point. Thus the prefactor 03(1, a(k)) has to be 
matched at tree- level to Eq. ([[i]), and the only solution to the leading logarithm is 



c 3 (l,a(fc)) =3(g R (k)Y 



l + O 



1 



(27) 



>g(/c/A QC D) 

Identifying the non-power-corrected term of Eq. ( |l8| ) to the purely perturbative coupling con 



stant, ci ( jp, a(fi) ) is known to be 



C\ ( ^2,Qt(/x) 



~ 5Ti,,pert(& 2 

Thus, we take from ref. || the following leading order results: 



1/2 



(28) 



'k 2 



7Q 

a(fc)\ 0o 



C2 ( -o,q;(m) 
\^ 2 



3^^ 



a(A:) \ 35 



a(/x) 



(29) 



where, as usual, 70 is defined by 



d In /i* 



lnZ(^ 2 



a(/x) 



(30) 



Eq. (p~8|) can be then applied to obtain a 



MOM 



g^/(47r) which, after the appropriate Wick 



rotation takes in Euclidean space the following form: 



a pcT t(k 2 ) 



1 + 



k 2 



In 



A 



2Q±2Q_i 

00 



/3 



with 



A) ivg - 1 



70+70 



(^4 2 )ij iM being now the Euclidean condensate; and 



(31) 



(32) 



«o = 7o - 7o , 

as immediately follows from Eq. (^) if kq is defined from 



(33) 



d 



lnZ K ( M 2 ) 



, „_,. w . , - -Ko I ... 

cZln/x 2 47r 
The following perturbative coefficients for the flavourless case, 



(34) 



A) = 11, A = 51, 70 = y 



(35) 



are universal, while coefficients in the MOM scheme, fa was computed in ||, and 71 and 72 in 




29 

ft ~ 4824. , 71 = — , 72 = 960 . (36) 
o 

In a recent pap er || we computed 70, from diagrams identical to these in fig. [2] (a) and (b); kq, 
defined in Eq. (|33"|) can be computed from diagrams (c) and (d) in fig. ^ to give: 

„ 3N C 9N C f . 

70 = — K0 = -l36-- (37) 

We have proceeded to evaluate the diagrams for kq in total analogy to the calculation of those 
for 70. For instance, we turned the diagrams to be infrared-safe by considering a momentum flow 
incoming to the local operator. The details of the procedure can be found in ref. ||. It should 
be noticed that kq << 70 and that, in practice, kq/Po — works as a good approximation to 
simplify Eq. (|3l]). In other words, the scheme given by vacuum insertion factorization reveals 
itself to be coherent: leading logarithm corrections violating the factorization induce a very small 
running with renormalization scale for the factorized tree-level Wilson coefficient. 

An important point is nevertheless to prove the assumption to work in order to obtain a 
good estimate of the Wilson coefficient for the asymmetric three-gluon Green function. To dig 
into this question, we have performed the same combined fits shown in ref. || for two and three 
points Green functions, at three loops for leading Wilson coefficients, to match lattice data taken 
from refs. fg, 13]. The Euclidean OPE formula for the two point Green function is, from ref. ||, 



7 MOM/l2 _\ _ 7 MOM/„2 „\ „ I k 2 , x 

7Q+7Q 

k \ Pq 



-1 



1 + RQi) (In-) T2 I ; (38) 



where 



«t M (M 2 ,a) 
and 



7MOM/L2 _\ 

^Latt ^ >«J _ fe 2 G g) (A .2 jM 2) + (a 2 ) , (39) 



7Q +7Q 



= AW^I) K) - (40) 

The coefficient cq(j^ , a(/j,)) is taken now to be expanded at three loops in terms of the MOM 
scheme a(k); i.e. it will verify the same differential equation as 7(a) in Eq. (p6|) with the 
boundary condition which is apparent from Eq. (|29|) . In the following, all the scale-dependent 
quantities will be shown at \x = 10 GeV. Furthermore, we have checked that both the ratio of 
gluon condensate estimates and A^jg- indeed do not depend on this last momentum scale. The 
quality of the fits as a function of a free exponent of In (-v) in Eq. (|3l|) has been explored and 
the results are shown in fig. |3|(a). We can conclude from these results that the approach given 
by vacuum insertion factorization provides a good estimate of this exponent. 
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Figure 3: (a) Quality of the fit as a function of the exponent r = 1 — 7 "^ 7 ° in Eq. 3J_. The dot stands 



for the value of r computed in this paper. In figure (b) is shown the fit of our Lattice data for oi^j^ [ at 
three loops with the calculated anomalous dimension. 





This work 


symmetric 3-point 


^MS 


260(18) MeV 


233(28) MeV 


{ > prop 


1.39(14) GeV 


1.55(17) GeV 


y^*A , „ 

<- J alpha 


2.3(6) GeV 


1.9(3)GeV 



Table 1: Comparison between results obtained for the three loops fit in the present work and in a 
previous one || 



The results for two particular values of the exponent r in Fig. ||(a) are of interest for the sake 
of comparison. The case r = — 1 corresponds to the formula proposed in ref. to be matched 
to the lattice data: a perturbative three- loop formula + a term c/p 2 , c being a constant. On the 
other hand, had we neglected the leading logarithm contributions, 7o = 7o = 0) the exponent 
would be r = 1. It can be seen from Fig. ||] that both values of r generate rather less good fits 
to the lattice data. 

Then, Eqs. can be used to perform fits at two and three loops for the leading 

Wilson coefficients in order to estimate, from asymmetric three-gluon Green function, the gluon 
condensate. The results of such a fits, plotted in fig. ||(b), are: 



- ( -f^ = 3.65(4) Ajrjs = 283(15)MeV X 2 = 1-95 (41) 



1 prop 

for the two loops fit, and 



1 n 



{ ^ i^*} alpha 



alpha _ A _ = 260(18)MeV = 1-18 (42) 



prop 

for the three loops one. 

The impressive improvement from two to three loops suggests that the approach presented 
in this work permits a reasonable approximation to the Wilson coefficient. The ratio decreases 
to almost two ex's away from 1 and both estimates of A^g and of the gluon condensate turn to 
be close of the previous estimates obtained from the symmetric three-gluon Green function in 
ref . H (see table [l]) . The scheme in this work and that of ref. || differ only by the kinematics 
of the renormalization point. Such a different renormalization for Green functions implies a 
different renormalization of the gluon condensate. However, the discrepancy for estimates of 
{A 2 )r in the two works is expected not to be important^, as indeed can be seen in table |l|. 
The comparison of these estimates, mainly the ones obtained from gluon propagator, strongly 
supports the claimed rather large contribution from the A 2 condensate that might be in 



connection with the tachyonic gluon mass scale studied in ref. 

A negative hint regarding to previous results from the symmetric three-point Green function 
is nevertheless the higher central value of the ratio in Eq. (42) (1.2 in ref. ||). In principle 



two possible sources of discrepancies could be expected: either three loops is still insufficiently 
accurate for the estimate of the perturbative part in the MOM renormalization scheme, or there 
is a deviation from the assumed vacuum insertion (or factorization) approximation. Both effects 

would have a direct impact on the bigger ratio we obtain for asymmetric MOM scheme. The 
very good agreement between the gluon condensates estimated from gluon propagator previously 
discussed seems to point out the factorization breaking as the major contributing factor. Still the 
ratio in Eq. (|42| ) is only about two sigmas from 1, which is in our opinion a rather encouraging 
result. The (A 2 ) deduced from the propagator is in fair agreement with previous estimates, 
even though it is biased by the factorization hypothesis through the fit of A^jg- which combines 
the propagator and as- This good result of the propagator as well as fig. ||(a), suggests that 
our formula describing the power corrections to as up to the leading logarithm yields a good 
approximation of the exact one. 

A two-sided goal is thus achieved: 

i) the results of ref. || turn out to be confirmed by the use of a slightly different renormal- 
ization scheme. 

ii) The vacuum insertion factorization applied to condensates playing the game of the OPE 
for the asymmetric three-gluon Green function results in a compact prediction for its OPE power 
corrections. The coefficient of the power correction has been computed to the leading logarithm, 
and thus a most important source of systematic uncertainty for the estimate of A^jg in ref. |7| 
is eliminated. The latter is a positive feature because the (3 function is perturbatively known at 



four loops in the asymmetric MOM [14] and lattice evaluations, on the other hand, turn out to 
be statistically more precise in this last renormalization scheme. 

A last consequence of this work and those from refs. |l,Ph they althogether lead to conclude 
that the Green functions methods, three-gluon vertex in particular, provide us with a reliable 
and precise enough estimate for the running coupling constant and A^jg, once power corrections 
are properly taken in consideration. 



3 To estimate the discrepancy for the non-perturbative estimates of (A 2 )n we need to compute beyond leading 
logarithm corrections that is out of the scope of this work. 
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